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Abstract 



It has been shown by several authors that a certain class of composite operators with many fields 
and gradients endangers the stability of nontrivial fixed points in 2 + e expansions for various 
models. This problem is so far unresolved. We investigate it in the A-vector model in an 
1/A-expansion. By establishing an asymptotic naive addition law for anomalous dimensions we 
demonstrate that the first orders in the 2 + e expansion can lead to erroneous interpretations for 
high-gradient operators. While this makes us cautious against over-interpreting such expansions 
(either 2 + e or 1/N), the stability problem in the A-vector model persists also in first order in 
1/N below three dimensions. 
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1 Introduction 



Problems with high-gradient operators were first observed by Kravtsov, Lerner and Yudson in 
the Q-matrix model in a 2 + e expansion |l|, ||] . They noticed that a certain class of canonically 
irrelevant composite operators with 2s gradients acquires anomalous dimensions in first order in e 
that grow with s 2 and that finally make an infinite number of these operators relevant for s — > oo 
and any fixed finite e > 0. Wegner et al. have shown that this behaviour occurs also in other 
2 + e expansions in one-loop order: In the A^-vector model ||], the unitary matrix model [Q] and 
the orthogonal matrix model [|| . In fact this behaviour has been observed in all 2 + e expansions 
that have so far been investigated with respect to this question. Judging from these results one 
could argue that the conventional nontrivial fixed points of these models become unstable against 
an infinite number of high-gradient perturbations. This would have important consequences for 
the present understanding of many problems that rely on 2 + e expansions. 

A natural starting point to investigate the role of high-gradient operators in more detail is 
the universality class of the A"-vector model. Here various expansions are available: 2 + e, 1/N 
and 4 — e expansions. In Ref. || Wegner found a composite operator with 2s gradients and the 
full scaling dimension 

y = 2(l-s)+e(l+ S 4^)+0(e 2 ) (1.1) 



N 

in the 2 + e expansion of the nonlinear c-model 

S = 2Y / ddx {d ^ )2 ■ (L2) 

Here 7? = (tt%, . .. ,7Tjv) is the normalized field tt 2 = 1 of the nonlinear a-model. Operators with 
y > 0, y = and y < are relevant, marginal and irrelevant, resp. Castilla and Chakravarty 
have extended the calculation to two-loop order || 0] and obtained the term in order e 2 

V - 2d " •) + « (l + ^) + j0sp (§ + 0(1/.)) + <>(«■). (1-3) 

This two-loop result makes the stability problem even worse. For s 3> N/e Eq. (|L^) amounts 
to an infinite number of relevant operators. On the other hand it could be proven in Ref. || 
that no similar problem occurs in the 4 — e expansion of (c/> 2 ) 2 -theory. In the 4 — e expansion all 
anomalous dimensions in order e can only make composite operators more irrelevant, but never 
go in the "dangerous" relevant direction. 

Therefore somehow the behaviour seems to bend over between d = 2 and d = 4 dimensions. 
The obvious tool to investigate this question is the l/A^-expansion for 2 < d < 4. This is the 
programme of this paper. Along its lines two other interesting aspects of the A^-vector model 
emerge. First of all it is shown how the complete spectrum of anomalous dimensions for a 
subalgebra of composite operators in order 1/N is encoded as a straightforward diagonalization 
problem. This formally solves part of a programme by Lang and Riihl Q. Then an asymptotic 
naive addition law for anomalous dimensions is derived that has also been found in the 4 — 



e expansion [10, 11]. This leads to a "clustering" of anomalous dimensions. However, this 
property is lost in a 2 + e expansion. It will become apparent that this serves as a warning 
against using the first orders of the 2 + e expansion to make statements about the high-gradient 
limit s — > oo. 

Still unfortunately this is not the end of the story. From an older result by Vasil'ev and 



Stepanenko [12] we will see that the high-gradient problem persists for 2 < d < 3 dimensions 
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also after adding up all order in e in first order in 1/N. From their work one can see that 
powers (7 s of the auxiliary field a in the 1/iV-expansion interpolate between the high-gradient 
operators in 2 + e and composite operators (4> 2 ) s in the 4 — e expansion. For some reason this 
simple fact has so far remained unnoticed and we also want to use this paper to clarify this point. 

Our paper is organized as follows. In Sect. 2 we discuss the renormalization of a certain class 
of composite operators to first order in the 1/iV-expansion. Eigenoperators and their critical 
exponents are encoded as a straightforward diagonalization problem similar to Ref. ||. In Sect. 3 
an asymptotic naive addition law of anomalous dimensions valid for a large number of gradients 
is found that has already been observed in the 4 — e expansion [10, 11]. However, it is shown 
that this structure remains unobserved in the 2 + e expansion which can lead to misleading 
interpretations. In Sect. 4 it is shown that results in first order in 1/N obtained by Vasil'ev and 
Stepanenko in Ref. [jO] allow the interpolation of the high-gradient behaviour between two and 
four dimensions. Sect. 5 contains the conclusions and suggestions for future work to solve the 
stability problem. 



2 The spectrum of anomalous dimensions in order 1/N 

2.1 The operator subalgebra C sym 



In this paper we shall mainly deal with composite operators Vl/^ 1 ^ 2 i" n (sc) constructed from n 
fields (f>A (A = 1, N) and their derivatives which are symmetric and traceless both with respect 
to their spatial 0(d) and internal O(N) indices. We denote this important class of composite 
operators by C sym , compare also Ref. fjD|. 

It is convenient to introduce "scalar" composite operators 

* (x, u, t) = *t:£ n {x)u h . . . uHa, ...t An , (2.1) 

where u and t are <i-dimensional and iV-dimensional vectors, resp. The original composite 
operators ^f^lf 2 '^" (x) can be reconstructed from ty(x,u, t) by differentiation with respect to u 
and t. 

Now ^Sf(x,u, t) can generally be expressed as 

(2.2) 



d d 

$f(x,u,t) = tp(-—,..., - — )$(x;u;t;ai,...,a n 
oai oa n 



a,i=0 



where 



<S>(x;u;t; ai ,...,a n ) = P{u,d v )T(t,d s ) ]J s Am <j> Am (x + a m v) 



m=l 



(2.3) 

v,s=0 



Here P and T are projectors on the subspace of traceless and symmetric tensors. The action of 
P (resp. T) on the operator ty(x,u, t) amounts to the elimination of all terms proportional to 
u 2 (resp.t 2 ). The explicit form of the projector P,T can be found in Ref. [14], but will not be 
used in the following. 

In this manner Eq. (|2.2j ) establishes a one to one correspondence between a composite operator 
ty(x, u, t) G C sym and a symmetric polynomial tp(zi, . . . , z n ), which we call the coefficient function 
of the operator ^(x, u, t). 

Now we briefly discuss the important subclass of conformal operators in C sym that will play 
a central role in our further investigations. The operator '^^...^(x) (for a moment we omit 
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the internal indices) is called conformal if it transforms as a tensor under conformal coordinate 
transformations, i.e. 

A 1 

<*a*ii,...,i,(aO = (ad + ^da)^ iu ... A (x) + Y / d ik a s ^i u .. ! s^ il (x) (2.4) 

d k=i 

at 

x — ► x'(x) = x + a(x), 5 a 0(x) = (a(x)3 + —da(x))(f)(x), 
with Qj(x) = tii + Wj fc a; fc + Axj + (x 2 5 ik - 2xiX k )c k . 

Here A and A$ are the conformal dimensions of the elementary field 4>(x) and the operator 
^h,...,ii( x )) r esp. Note that for conformal operators from C sym one simply has A^ = (scale 
dimension) - (number of spatial indices).) 

A straightforward calculation shows that the coefficient function tp of conformal operators 
must satisfy the following equation 

n 

(*mC + Ad Zr M( Zl , ...,z n ) = 0. (2.5) 

m=l 

To find all solutions of this equation let us define the one to one map ip — > tp: 

^ 1 poo n 

^(z 1 ,---,z n ) = Tip(z l ,---,z n ) = — r^r / Hdt m t*- 1 e- tm ip(t 1 z u ...,t n z n ). (2.6) 

\ ' m=l 

For the function ip Eq. ( |2.5| ) reads 

n 

Y,d Zr J( Zl ,...,z n ) = 0, (2.7) 

m=l 

that is ip must be translationally invariant. Thus the space of conformal operators with n fields 
and I derivatives is isomorphic to the vector space C(n, I) of symmetric, homogeneous of degree I 
and translationally invariant polynomials of n variables. The dimension of C(n,l) can be easily 
calculated from the generating function D(n,x) 

oo -. 

D(n,x) = ]TdimC(n,l) x 1 = — - — -. (2.8) 

1=0 llp=2U- xP J 

2.2 Calculation of Feynman diagrams 

To describe the renormalization of operators from C sym in the nonlinear <r-model we follow the 
scheme developed in the papers [|T^, 15]. The regularized action for the theory under consideration 
has the form 

S = --d^d^A - \aK u a + X -M"o<\> A <\> A + ^M^aKa . (2.9) 

Here uj and M are regularization parameter and regularization mass, resp. The kernel K is 
determined from the requirement of the cancellation of self-energy diagram contributions of the 
sigma field (see Fig. 1); K(x) = —N/2G 2 Ax), where G^{x) is the propagator of the 4> A field 
(<<P A (x)^ B (y)>=S AB G 4> ( X ))^. 
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The regularized kernel K u is defined as K ul (x) = K(x)x 2uj . 
The explicit expressions for the propagators and G a = K~ x in momentum space are 

G^p) = p- 2 , G (p) = B p- 2 ^- 2+ "\ (2.10) 

where fi = d/2 and B = -2 • 4^ +w vr^r(2/i - 2 + u)/ (NT 2 (fi - l)r(2 - fi - w)). 

In the calculation of Feynman diagrams the divergencies appear as poles in uj which therefore 
corresponds to e in the dimensional regularization scheme. It should be noted that the nonlinear 
cr-model is not multiplicatively renormalizable 15]. In general its critical exponents cannot be 
expressed via the renormalization constants Zi but must be determined (employing conformal 
invariance of the theory [15]) directly from the renormalized Green functions by looking at the 
long distance logarithmic contributions. However, in Ref. [12] it was shown that in order to find 
the critical exponents in first order in 1/N it is sufficient to calculate the divergent pole term of 
the Feynman diagrams. A detailed discussion of this point can be found in fi"2]| . 

Thus to find the critical exponents in first order in 1/N we must extract the divergent part 
in the corresponding diagrams. 

For the operators from C sym only the diagram in Fig. 2 contributes in first order in 1/N. 
Diagrams containing cycles of <j> lines are zero due to the O(N) tracelessness of the operators 
in question. Diagrams describing transitions didk — > <5ifccr vanish due to the tracelessness in the 
spatial indices. 

As usual one has to consider the renormalization of the set of mixing operators. Renormalized 
operators are defined as 

The critical exponents are obtained from the mixing matrix Q via 

xt{g) = n^5 rr> +2Q rr \ (2.12) 

where n is the number of fields cj) in the operator \I'r( < / ) ) and 7^ = rj/2 is the anomalous dimension 
of the elementary field 4> with 

4 T(2^ - 2) 2-fj, 2 

71 = n r^-i)r ( 2-^) foTPi) + 0{1/N ] ■ (2 - 13) 

It is convenient to fulfill the renormalization procedure in terms of the generating func- 
tion (|2.3|). The analytic expression corresponding to the diagram shown in Fig. 2 is 



r ]2fii i(a.i-a k )(u,l) 

i(u,piai+p k a k ) / " ^ (o -1 A) 

J (2tt) 2 M p(Mr-2+w)(l +p .)2( pk _ l y> ^ - > 

where (•,•) denotes the scalar product. The calculation of this integral does not cause any 
problems and we give the final answer for the renormalization of the generating function at once 

[$(x;u;t;ax, ...,a n )] R = $(s;u;t;ai, ...,a n ) + ^-y^H ik ^(x;u;t]ax, ...,a n ) , (2.15) 

UJl\ . , 

i<k 

where 

, x T(2u-2) 
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and 



H ik $(..ai, .., a k , ••) = l/2(h ik + h ik )$(..a,i, .., a k , ..) = 

l / 2 ( I IT da m <5( V a n - l)a%~ 3 ®(.., aid* + aiojt, ..,a 3 afc + a 3 ai, ..) + (i <-> fc)") (2.17) 
^° m=l n=l " J 

with the notation a m = (l—a m ). On the rhs of this integral equation the i— th and fc— th arguments 
are replaced, the other arguments (denoted by dots) remain unchanged. The first term in the 
brackets in Eq. ( |2.17| ) det ermines the action of the operator h^ k on the function $(.., o>%^ .., d k ^ • •), 
while the representation for hi k is obtained from the latter by interchanging o« <-> a k . Although 
on symmetric functions hi k and h ik are identical, it will be useful to distinguish them in the 
sequel. 

To avoid misleadings we stress that ( 2.15| ) should be understood as an equation for generating 



functions, i.e. as an equation for the coefficients in the power series expansion in of the lhs and 
rhs expressions. Eq. ( 2.1 7^ ) also requires some comments. Due to the factor o^ -3 the integrand 



has a nonintegrable singularity and the integral in rhs of Eq. (2.17) must be understood as the 
analytical continuation from the region [i > 2. As a matter of fact a^ -3 * s always multiplied by 
some polynomial in and the integral can always be expressed in terms of T-functions. 

2.3 Diagonalization problem 

To obtain the spectrum of critical exponents we must diagonalize the mixing matrix Q, i.e. 
find the multiplicatively renormalized eigenoperators. Taking into account Eqs. (2.2), ( p. 11 ), 



( 2.15| ) it is easy to derive the eigenvalue equation for Q. Indeed, using the Fourier representation 



for (j)(x + diu) (4>(x) = J d 2>1 q 4>(q)e tqx ) one immediately finds that the operator ^(x,u, t) is 
multiplicatively renormalized if its coefficient function satisfies the following equation 

Xip(zi, ...,z n ) = ^Hikipiz!, ...,z n ) = l/2^2(h ik + h ik )ip(zi, . . .,z n ), (2.18) 

i<k i<k 

where hi k again has a two-particle form 

f 1 3 3 3 

h ik ip(..Zi..z k ..) = / Y[ da m 5(^2 a n - l)a%~ a\Zi + a 3 z k , .., axZi + a 3 z k , ..) . (2.19) 

^° m=l n=l 

Again on the rhs of this equation the i— th and k-ih arguments are replaced in the manner shown. 
The expression for hi k is obtained from ( 2,19| ) by interchanging Zi <-> z k . The critical exponent 
of the corresponding eigenoperator is 

x = n{n - 1) + I + n| + + 0(1/A^ 2 ) (2.20) 

corresponding to a full scaling dimension y = d — x. Here and in the sequel n is the number 
of fields <f) and / is the number of derivatives in the eigenoperator. The function v((J,) has been 
defined in ( 2.16 ) and rj was given in ( |2.13| ). 



Eq. ( |2.1&j ) has a number of nice properties. But the eigenvalue problem is more manageable 



if written in terms of the ifj(zi..z n ) functions. Under the transformation (2^) the lhs and rhs 
of ( |2.18 ) transform into 



AV>(£i, ...,z n ) = ^flife^i, ... ,z n ) = H^{z\, ...,z n ) (2.21) 

i<k 
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with the operator defined in Eq. ( 2.17| ). Some technical details of this calculation are given 
in Appendix A. 

Now it is easy to show that the symmetry group of H is SL{2, C). We define the action of 
the SL(2, C) group on functions ^p(z±, . . . , z n ) in the following manner 

S(g)^( Zl , . . . , z n ) = n(cz ; + d)" A # ?*±±) (2.22) 

Vczi + d cz n + d/ 

where g £ SL(2, C), g = ^ ° ^ ^ and ad — be = 1. It is easy to check that H and S(g) commute 

for all g with A = fi — 1 in our case. The group SL(2) has three generators S , S+ , S_: 

n n n 

S = 5>i<9 Zi + A/2), S_ = -5> ZiJ S+ = ^(zf9 Zi + Azi). (2.23) 

i=l i=l i=l 

All these generators commute with H and obey the commutation relations 

[S, S_] = -S_, [S, S+] = +S+, [S+, S_] = 2S. (2.24) 



Using Eqs. ( 2.24 ) it is straightforward to see that the eigenvalue problem for H can be restricted 



to the subspace of functions annihilated by S_: 

Sf n = {fe, • • • , *»)|S_# = 0; = (1 + nA/2)V4 • (2.25) 



According to (|2.7| ) the operators in S_ are just the conformal operators transforming like 
Like in Ref. || we therefore observe that the spectrum of critical exponents is generated by the 
conformal operators. All other eigenfunctions of H are obtained by application of S+ on the 
eigenfunctions in , (I, n = 0, 1, . . .). Thus we have to solve the eigenvalues problem (|2.21|) for 
the operator H on the space of functions ip(zi, z n ) that are symmetric, translation invariant 
and homogeneous polynomials of degree I in n variables. 

Another very useful property of H is its hermiticity with respect to a suitable scalar product 
defined on the space of homogeneous polynomials as 

<$i\$2 >=M9z 1 ,...,d Zn )Mzi,...,z n )\ Zt=Q . (2.26) 
Remember that ip and if) are related via transformation (|2.6| ). To prove the hermicity of H with 



respect to the scalar product (2.2 



< ifalHfo >=< Hfalfa > , (2.27) 



it is sufficient to note that Hif.tp(zi, . . . , z n ) —> Hikip(z\, . . . , z n ) under the transformation (|2T 

An important immediate consequence of this hermicity is the reality of the spectrum of critical 
exponents, i.e. all anomalous dimensions of eigenoperators in C sym are real numbers at least in 
first order in 1/N. 

2.4 Double expansions 



A closed analytical solution of Eq. ( 2.21 ) seems not possible except for the trivial case n = 2 where 
the form of the eigenfunctions is already fixed by conditions imposed on the functions ip(zi, Z2) 
alone. Indeed there is only one symmetric translation invariant polynomial of degree I = 2L: 



4>i{zi, Z2) = (z\ — Zq) . The corresponding eigenvalue of ( 2.21 ) is A; = !/((/ + // — 2){l + [i — 1)) 



leading to critical exponents ( |2.20| ) in agreement with known results [1 
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There are two important cases when the situation is simplified and the problem becomes 
more tractable. Let us consider the eigenvalue problem in the limiting cases where d = 2 + e or 
d = 4— e, in other words we want to hold the leading terms in e in the anomalous dimensions. This 
corresponds to a double expansion in 1/N and e keeping the leading nontrivial terms in Eq. fl2.20p . 
The multiplier v(fi) is in both cases proportional to e 2 , whereas A determined from Eq. ( |2.21| ) has 
a simple pole in e. Extracting the pole contribution from the expression Hip(zi, . . . , z n ) (defined 
by Eq. ( |2.17 )) one obtains the following formulae for the critical exponents with 0(e 2 ) accuracy 
in d = 2 + e 



X = n 2V + N] +l+ 2N X2+e + ° {e) + °W N ) ( 2 - 28 ) 
and d = 4 — e dimensions 



x = n f 1 - -\ + I + -A 4 - e + 0(e 2 ) + 0(1/N 2 ) , (2.29) 
where A2+ e and \i- e are given by the eigenvalue problems 

\ 2+e $(zi, ...,z n ) = ($(.., Zi,..., z h ..) + $(..,z k , . . . ,z k , ..)J (2.30) 
i<k 

and 

A 4 _ e ^(zi, . . . ,z n ) = 2 V / da>Tp(..,Zia + (1 - a)z k , . . . + (1 - a)z fc , ..) . (2.31) 
i<k Jo 

Again the i-th and A:-th arguments on the rhs of these equations are replaced in the specified 
manner. 

These double expansions can be compared with known results, which serves as a consistency 
check of our calculations. Formula ( 2.3l|) has already been obtained in Ref. [13] in a 4 — e ex- 
pansion. In fact, one can also check in a lengthy calculation that the term in order e 2 /N from 
Ref. [10] is in agreement with our 1/iV-results. 

In the 2+e expansion the complete spectrum of critical exponents in order e has been obtained 
by Wegner ||. The solution of the eigenvalue problem ( p. 30 ) is straightforward and yields 

A2+ e = m(2n — m — 1), m = . . . n . (2.32) 

Here m is the maximum number of elementary fields without any gradients acting on them in 
any term of the respective eigenoperator. This is in agreement with the results of Wegner ||] 
restricted to our subalgebra C sym as can be checked easily. 



3 An asymptotic naive addition law for anomalous dimensions 
3.1 Outline of the proof 

Since we have the exact solution of the eigenvalue problem for operators from C sym in first order in 
the 2 + e expansion, it is natural to try to understand how large the corrections from higher orders 
in e are. Indeed, the calculation of critical exponents of high-gradients operators in C sym will 
lead to the observation that an appropriate expansion parameter is not e but some combination 
of e and Later we will see that a suitable combination is e.g. e • InZ. 

The 1/iV-expansion, in which the higher orders in e are summed up, is very suitable for this 
problem. However, due to very fast growth of the size of the mixing matrix with increasing I it 
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seems hardly possible to obtain the analytical solution of the eigenvalue problem Eq. ( |2.21| ). But 
the following interesting quantitative prediction for critical exponents of high-gradient operators 
can be done. 

Theorem: Let be the spectrum of H = J2i<k^ik restricted to the subspace S™'\ i.e. 
the subspace of symmetric homogeneous polynomials of degree I in n variables. Then for any 
two eigenvalues Ai G A^ 1 and A2 € A^ 2 2 corresponding to different conformal eigenoperators 
the following statement holds: There is a number L such that for all[] / > L one can find an 
eigenvalue A; 6 A* +n such that 

|A,-Ai-A 2 | <C/l A . (3.1) 

Here C is a fixed constant and A = fj, — 1. 

Hence any two anomalous dimensions for composite operators with n\ and ri2 elementary 
fields generate a limit point given by the sum of these anomalous dimensions in the spectrum of 
composite operators with n = n\ + ri2 fields. The approach to this limit point is given by (|3.1| ). 
Therefore Eq. ([O]) is an "asymptotic naive addition law" for anomalous dimensions. The limit 
point structure in the spectrum leads to a clustering of anomalous dimensions. For more details 



we refer the reader to Refs. [10, 17 1 where this property was discussed in the framework of the 
4 — e expansion. 



The general scheme of the proof of the theorem is the same as in the paper [17]. So here we 
will only discuss the main ideas and some technical modifications as compared to [17]. The full 
proof can then be constructed along the same lines as in Ref. [17]. 



To prove the unequality (|3.l| ) it is sufficient to show that in the subspace S 1 ^ 1 there exists a 
vector ^5(2:1, ■ • • , z n ) such that 

[|(if-A 1 -A 2 )§ s ||/||*5||<C/Z A . (3.2) 

Here, as usual, ||$s|| 2 =< $s|$s >• 

We begin with the construction of the test vector tyg. If tpi(zi, • • • , z ni ) and ip2 (zi , • • • , z n2 ) 
are the eigenfunctions of H corresponding to the eigenvalues Ai and A2, we define ^$ m the 
following way 

$5(21, • • • , zn) = T* s (2i, • • • , z n ) = Sym T#(zi, ■ ■ ■ , z n ) , (3.3) 

{zi,...,Z n } 

where Sym denotes symmetrization with respect to z%, ■ ■ ■ , z n . The function ^(z\, ■ ■ ■ , z n ) reads 
*(x,y) = °k( x i + ■■■ + x ni ) k (yi + ■■■ + y n2 ) l ~ k Mx)My) , (3-4) 

k=0 

where c k = (-l) l - k C l k r(S 1 )r(S2)/r(k + S 1 )r(l-k + S 2 ), Si = mA + 2Ri and C l k is the binomial 
coefficient. Also, for brevity, we use the shorthand notation ipi^i^) = ^1,2(^1, ""> x ni 2 ) and 

^(x,y) = • • ■ ,x ni ,yi, ■ ■ -,y m )' 

A straightforward calculation shows that ^f(z) defined by Eq. (3.4) satisfies Eq. (|2.5|). Con- 



sequently, the functions ^(z), ^s(z) are translation invariant. The degree of the polynomial 
^s(z) is obviously equal to / + R\ + Bq. 

1 If Ai and A2 correspond to the same eigenoperator, this statement holds only for even I > L. 
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To prove the unequality (^) for $5 given by Eq. ( ^^ ) one needs to calculate \\$s\\ and 
\\(H — Ai — A2)^ , 5||- To do this we write the functions *$>s an d i n the following more suitable 
form 

m n 2 

(3.5) 



(x, y) = i)i(x)i) 2 (y)K{d a , 8 b ) exp{a(J2 x i) + Vj)} 

i=i 3=1 

where if (a, b) = J2k=o Cka h b l ~ k . 



a,b=0 



Respectively, applying the transformation (2.6) to Eq. ( |3.4| ) yields after some algebra 
$(x,y)=/ dad(3 [ DuDvf(a,n 1 ,R 1 )f(P,n 2 ,R 2 )[a(ux)-(3(vy)] l ip 1 (ux)tlj 2 (vy), (3.6) 



where f(a,n,R) = (1 - a) 3 -^ 1 ^- 1 ■T{S)/r(R)T n (A); Du = U?=i <foi«f ~ 1<y (Di = 
J2iUiXi &nd*p(ux) = ip(uixi, . . . ,u ni x ni ^. 

Now the calculation of the norm of ^5 can be done in the following manner: 
1. Using formulae ( [2.26 ), (3.3) and taking into account that the functions ipx i2 (z) are fully 
symmetric one obtains 



"2 



= (Cn 1+ n 2) -l £ CJI^W , (3. 7) 
k=0 

where A, (fc) = *(9* s • • ■ j 2/fc> ^fc+l) ■ ■ ■ > x ni) x li ■ ■ ■ i x k: Vk+li ■ ■ ■ 1 2/712)- Without loss of gen- 

erality we now assume n 2 < n\. 

2. Using formulae ( |3.5[) , (|3.6| ) and the translation invariance of $ one can obtain the following 
expression for A^ 

Af ] = N(l) f 1 dad/3 ^ DuDvf(a, rn,Ri)f(fi, n 2 , R 2 )A l {a, f3, u, v)T{a, 0, u, v) , (3.8) 
Jo Jo 

where A{a, f3, u, v) = [a{U k - U k ) + f3(V k - V k )}/2; U k = J2 k =Q «i! V* = E*=o «»; 
N(l) = K(d a ,d b )(a - b) l \ a=b=0 and 

R1+R2 

J r (a,..)= ^2 C\A~ l d] l ^i{d x —a(3v 1 d x —aau)il}2{dy+aau 1 dy—al3v)^i{uy,ux)il} 2 {vx,vy) x=1 / 2 ■ 
i=o y=-i/2 

Here ipi(x,y) means ipi(xi, ..,x k ,y k+ i, "jj/nj and we mention again that x = (1 — x). 

3. All the dependence on / in expression ( |3.8j ) for -A; , except the trivial one contained in the 
multiplier N(l) (we do not need its explicit form), is due to the factor A 1 in the integral. It is 
evident that the main contributions in the integral come from the regions in which |«4(...)| — 1: 

{A ~ 1 a ~ p ~ [4 ^ U fc ~ 0; }; {.A =± -1 a ~ /? ~ C/ fc ~ V k ~ 0} . 

Then after long but straightforward calculations one finds that the leading contributions in the 
integral Q from these regions behave as cf ] /l 2kA for A ~ 1 and c5/l(" 1+ " 2 - a ) A for ,4 ~ -1. 
The explicit expressions for the constants are 



^ = rnlf^-if^) / / Pfe+i^] ^ 2 M, (3.9) 



^ = (- 1 ) | r( ^ ) ^ l) / (^(«,0)^(u,0)) / [23f«] (^,0)^,0)), (3.10) 
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where [D%u] = YYU ^ ^ <*(E2 «i - 1). 

"~ • (h) 

Thus one can see that the main contribution to the norm of is given by A) when k = 

(k) 

or 2k = n\ + n^. It can be shown that in these cases the exact answer for A\ can be obtained 

(k) 

and it coincides with the answer for the leading term of the asymptotic behaviour of A\ . Put in 
another way, the functions with arguments arranged in different ways are almost orthogonal 
to one another for large I. 

Eventually the expression for the norm of *$>s is 

||$ s || 2 = N{l){C^)- l \\M 2 \\MH^ + (-iy<W(! + 0(l/l A )) ■ (3.11) 



Our next task is to obtain an estimate from above for £ (I) = \\(H — Ai — A2) 1 9/,s||. Instead 
of trying to derive the exact asymptotic behaviour for this quantity, we prefer to give a rough 
estimate. This, however, reproduces the true dependence on I but in a much shorter calculation. 

1. First of all, taking into account that H = J2i<k ^ik is symmetric under transpositions of its 
arguments and using Schwartz's inequality one gets 

£(l) = \\(H - Ai - A 2 )*<j|| < \\(H - Ai - A 2 )$|| . (3.12) 
Remember that ^?(z) = ^(z)(z\ ■ ■ ■ z n . 

2. It can be easily shown that (z\ + ...z n )ip(z) — > S+^>(z) under transformation fl2.6| ). Then taking 
advantage of the explicit form of ^(z) (Eq. (J3.4|)) and the SL(2, C) invariance of H one obtains: 
Y^i<kHik^{z) = \\^(z) and J2n 2 <i<k Hik^(z) = \2^(z). Taking into account the symmetry 
properties of ^?(z) and again using Schwartz's inequality we find 

£{l) < nwz < §{x, y)H 2 ( Xl , yi )${x, y) > . (3.13) 



It is possible to simplify Eq. ( 3.13| ) even more. Let us remember that the eigenvalues of = 
l/2(h ik + h ik ) are A/ = [(l + fi- 2)(Z + fi- l)]" 1 I = 0,2 ,4... One sees that for 1 < \i < 2 the only 
negative eigenvalue is Ao- All the other eigenvalues lie in the interval (0, 1). Thus if Po is the 
projector onto the subspace of eigenfunctions of H(x, y) with A = Ao one has < ty\H(xi, yi)]^ >< 
< $|(H-A P )|$ > +Xl < $|P |$ >< < *|fc(a?i,yi)|* > +A (A - 1) < $|P |$ > 
The projector Pq can be written as 

P ^(x, y) = B-\A, A) f 1 ds(s(l - s))^ 1 ^ + (1 - s)y, sx + (1 - s)y) . (3.14) 

J o 

Here B(x,y) is Euler's beta function. 

3. The calculation of < ^ / [Po| x I / > and < 5 r |^.|^ r > is carried out in a similar manner, so we 
consider the latter only. First of all it is useful to rewrite the expression (|2.17| ) for in a 
well-defined form for 2 < d < 4 (compare Eq. ( 2.17] ) and the comments to it). 



hi k i){x,y) = (A - l) -1 / dsdtt(l - t)^ 1 (/i + td t )^((l - st)x + sty,tsx + (1 - st)y). (3.15) 

J o 

Furthermore, in full analogy with the derivation of Eq. ( p.8| ) we obtain 
Nil) f 1 

< ^\h(xi, y i)y& >= w / dadf3dsdtDuDv f(a,n ll R 1 )f(f3,n 2l R2) 
A — 1 Jo 

x t(l - i) A-1 (ji + td t ) A l (a, (3, u, v, t, s)Q(a, (3, u, v, t, s) , (3.16) 

where A = [a(l/2 — u^ts) + (3(1/2 — v\ts)]. The function Q is a polynomial of degree Pi + P 2 
and has a structure like T in Eq. (|3.8|) Again one concludes that the dominant contribution 
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to the integral comes from the region where \A\ — 1. At first sight the terms arising from tdt 
in ( p. 16 ) give the leading term of the asymptotic behaviour as / — > oo, since differentiation of 



A 1 produces additional power of I. But these additional powers of I come in pair with "small 
variables" ~ (1 — .4) that lead to the observation that the "^-term" and the "£<9i-term" in 



Eq. ( 3.16 ) have the same order in /. Both these terms can be treated on an equal footing. We 
consider the first of them. The condition 1 — A < e leads to a ~ (3 ~ u\ts ~ v\ts < e. Then 
in this region A (1 — a/2)(l — 0/2)(l — mts — vits)) + 0(e 2 ). Function Q in its turn can be 
bounded by Const ■ l R i+ R i, Using this the integration with respect to U2, -u ni , u 2 , --V n2 and a, (3 
can be performed easily giving a factor l~( R i+ R v . Thus the integral under consideration can be 
bounded from above by the following expression 



C(A)N(l) I dudvdsdttil-t^u^-^^v^-^-^l-ust-vts)) 1 





After tedious but straightforward calculations one shows that the leading term of the asymptotic 
expansion of this integral behaves as Z~ 2A . So the final answer for £ (I) is 

£(!) < CN(l)r 2A . (3.17) 
It is clear that the unequalities (|3.2|), ( |3.1| ) now follow immediately from Eqs. ( 3,11| ), ( |3.17| ). 



3.2 Failure of the limit point structure in the 2 + e expansion 

From the results in the previous subsection it is clear that the approach to the limit point close 
to two dimensions is like l~ t l 2 , where I is the total number of gradients. Hence the behaviour is 
non-uniform for d £ [2,4] and only uniform for d € [2 + 5, 4] with some fixed 8 > 0. A suitable 
expansion parameter for high-gradient operators close to two dimensions is therefore necessarily 
some combination of e and I, e.g. the most straightforward choice is e • In/. 

Let us make this more explicit by looking at the spectrum of conformal operators with 
n = 3 fields. The anomalous dimensions obtained by diagonalizing (2.21) are plotted in Fig. 3. 



Let us concentrate on the curves approaching the limit curve labeled (ft 2 <g> <j) given by the (trivial) 
sum of anomalous dimensions of <f) A and hA 1 A 2 ( t )Al( f >A2 ■ Here hA x A 2 is a symmetric and traceless 
tensor. First of all one can readily check that the approach to the limit curve as a function of I is 
m fact like l~ d / 2+l , with even I approaching from above and odd I from below. The approach is 
clearly non-uniform in the interval d E [2, 4]. Close to two dimensions the derivatives of the curves 
even diverge for I — > oo, corresponding to diverging terms in order e 2 /N in the 2 + e expansion. 

That the limit point structure cannot be observed in the 2 + e expansion is also apparent from 
the diagonalization of ( [2.30 ). We find the following possible eigenvalues for different numbers of 



fields n according to ( p. 32 ): 



n = 1 
n = 2 
n = 3 



A 2+e = 
A 2 +e = 0, 2 
A 2+e = 0,4,6 



This of course agrees with the values in Fig. 3 for d = 2. The naive sum of the eigenvalues 
belonging to 4> A and hA 1 A 2 ( t >Al( t )A2 with n = 1,2 is A 2 + e = 2, but this value is not contained in 
the above list of possible eigenvalues with n = 3 fields. Similar observations can be made quite 
generally in the spectrum of the 2 + e expansion, i.e. the limit point structure ( |3.1| ) remains 
invisible. 
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4 Vasil'ev and Stepanenko's results re— examined 



4.1 Interpolation of the high— gradient behaviour for 2 < d < 4 

In the operator subalgebra C sym discussed so far we have seen that one must be cautious if 
one uses e as an expansion parameter in d = 2 + e for high-gradient operators. This can lead 
to wrong conclusions about the limit of a large number of gradients I — > oo. Now, however, 
the problematic class of high-gradient operators giving rise to Eqs. ([!]), (|L3|) does not belong 
to C sym . The structure of the respective eigenoperators in the 2 + e expansion is || 

l 

lk(d+n ■ d^n) 2( - l - k \d + n ■ d + ^) k {d^ ■ d^) k (4.1) 

with some coefficients 7& (70 7^ 0) and d± = ^(d x T id y ). Hence these eigenoperators of the 
nonlinear cr-model ( |1.2| ) are O(N) and 0(d = 2)-scalar. 

In this section we use the opportunity to show that the operators ( |4.1| ) have in fact already 
been investigated in the first order of the 1 /TV-expansion for 2 < d < 4 by Vasil'ev and Stepa- 
nenko [12|. As this is not immediately apparent from their work, we want to clarify this point 
here. 



In Ref. 1 12] Vasil'ev and Stepanenko calculated the critical exponents of the composite op- 
erator a s made up of powers of the auxiliary field a from ( |2.9[) .p| They find the following full 
scaling dimension 

, (s - l)d(d - 3) - 2(d - 1) o N , . 

y = d-2s-i 1 s{d-l) y - \_ d " L + 0(N- 2 ) (4.2) 



with r\ from ( 2.13 ). If one expands this additionally in d = 2 + e one finds 

y = d - 2s + ^s(s - 1) + 0(e 2 ) + 0(V" 2 ) (4.3) 
or alternatively in d = 4 — e 

y = d-2s- 6-^s(s - 2) + 0(e 2 ) + 0(N~ 2 ). (4.4) 



One readily observes that (|4.3f) corresponds to the high-gradient operators (|4.1|) as the crit 



ical exponent is consistent with (1.1). In d = 4 — e Eq. fl4.4j) corresponds to the composite 
operators (4> 2 ) s with the full scaling dimension [18] 

y = d-2s-e^^ + 0(e 2 ). (4.5) 

As these critical exponents are unique in the respective expansions, this proves that powers of 
the auxiliary field a s interpolate between the high-gradient operators ( |4.l[) with 2s fields and 
2s gradients in d = 2 + e and the composite operators {4> 2 ) s with 2s fields and no gradients 
in d = 4 — e. This might appear surprising at first sight as both a s and (4> 2 ) s contain no 
gradients. However, the canonical dimension of <j> is x</, = d/2 — 1, which varies as 

d = 4 - e : x$ = 1 — e/2 
d = 2 + e : x <t> = e/2 



2 In Ref. ]l2| the auxiliary field was denoted by ip. 
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and hence allows to built gradients into the respective eigenoperators as d varies from four to 
two. Therefore the "most stable" eigenoperators ((p) s in d = 4 — e are linked with the unstable 
high-gradient operators ( |4,1| ), 

Eq. ( [4.2|) allows to answer the question where this behaviour bends over between two and 
four dimensions. Eq. ([4.2|) can be written in the form 



y = d - 2s - (s(s - 1) K(d) + O(s)) + 0(N- 2 ) (4.6) 
with 

m=W-*) Pm _$e% )TW) (■*-!)■ (4-7) 
K(d) is plotted in Fig. 4. It has the following general behaviour: 

K(d) > for d > 3 
if(d) =0 for d = 3 
K{d) < for d < 3 



As K(d) multiplies the quadratic term in s in fl4.6j ), this yields an infinite number of relevant 
operators in first order in 1/N for 2 < d < 3 in the limit s 3> N. Hence d = 3 is the dimension 
where this seemingly unstable behaviour changes into the stable behaviour of the 4 — e expan- 
sions ||. Of course this observation can be affected in any direction by the unknown terms in 
order 1/N 2 . 

4.2 Contact exponents 

As a short aside we briefly put the above observations into perspective of a scenario suggested 
by Duplantier and Ludwig |l9[| . They have argued that only field theories with negative contact 
exponents O < have the potential to describe multifractal scaling phenomena (leaving aside 
the use of analytical continuations like e.g. the replica method in field theory). A contact 



exponent [20] O = x k — X{ — Xj is defined from the short-distance OPE 



Oi(n) Oj(r 2 ) = |ri - r 2 \ Xk Xi Xj O k (r 2 ) + subdominant terms with x k > > x k . (4.8) 

Multifractal scaling requires [p| < 0, whereas in field theories one generally finds > @. 
E.g. in the 4 — e expansion one has for the symmetric traceless part of (f) Al . . . (f> An the unequality 
x n — nx\ = n(n — l)e/(N + 8) + 0(e 2 ) > in agreement with > 0. Similar for ((j) 2 ) s one finds 

x {$2)a ~sx f2 = 6s(s - l)^g + 0{e 2 ) > . (4.9) 



However, as we have seen in the previous subsection the behaviour (|4.9| ) changes qualitatively 
below three dimensions in first order in 1/N. For powers of the auxiliary field a s one deduces 



from Ref. |2| 



sx a = ^s(s-l)K(d) + 0{N- 



> for d > 3 

= for d = 3 (4.10) 
< for d < 3 
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with K(d) defined in (|4.7|) . Hence below three dimensions one finds the unusual behaviour for 
field theories 



< s x a for 2 < d < 3 



(4.11) 



in first order in 1/N allowing for negative contact exponents < 0. In the framework of the 
2 + e expansion this was already noticed in Ref. [19]. Notice that in Eq. ( 4.11| ) one is not forced 
to consider the limit s 3> N, (4.11) holds for arbitrary s and is in this sense a truly perturbative 
expansion in a possibly small parameter 4>. 

Obviously, all this makes sense only if the usual nontrivial fixed point does not become 
unstable for s 3> N, 2 < d < 3. Which cannot be answered conclusively for the time being as we 
have seen in the previous subsection. 



5 Conclusions 

In this paper we have investigated properties of high-gradient operators in the iV-vector model. 
This class of operators shows a seemingly unstable behaviour due to the large corrections from 
anomalous dimensions (1.1), (1-3) found in Refs. @,|6|, 0- Therefore it was the natural programme 



of our work to use an 1/iV-expansion for 2 < d < 4 dimensions to learn something about the 
behaviour of these operators as a function of d. Here we briefly sum up the main results. 

For the subalgebra C sym of composite operators introduced in Sect. 2.1 we have obtained the 
spectrum of critical exponents in first order in 1/N in a similar manner as in Refs. |||, [l3|, [Rj 
for 4 — e expansions: The spectrum is encoded as a straightforward diagonalization problem 
( 2.21| ) in the space of homogeneous, symmetric and translationally invariant polynomials. Let us 



mention that the spectrum of critical exponents in the 1/iV-expansion has been the subject of 
various papers by Lang and Riihl || that aimed at finding algebraic structures in d dimensional 
conformal field theories. In C sym our formalism gives the same results as theirs based on operators 
product expansions. 

In Sect. 3 we proved a similar limit point structure (or: naive asymptotic addition law for 
anomalous dimensions) as seen in Refs. |To| , |Tl|| in 4 — e expansions: Two eigenoperators with 
anomalous dimensions Ai,A2 and n\,ri2 elementary fields generate the limit point Ai + A2 in 
the spectrum of conformal operators with n\ + ri2 elementary fields. The approach to this limit 
point is like | 1-f V 2 where I is the total number of gradients. Such "limit curves" for anomalous 
dimensions (see Fig. 3) have already been observed in Ref. |T^] f| 

In the 2 + e expansion, however, this limit point structure cannot be seen as the approach 
is like r e / 2 if one uses e as the expansion parameter (for a detailed discussion see Sect. 3.2). 
Unfortunately this observation could not solve the stability problem of 2 + e expansions as 



the problematic class of operators leading to (|l.3|) does not belong to C sym . Still, besides the 



interesting non-perturbative structure also its failure in the 2+e expansion teaches the important 
lesson that e might not be the suitable expansion parameter for high-gradient operators. E.g. 
one would also be tempted to read of a stability problem in order e 2 in Fig. 3 close to d = 2 for 
curves with even I tending to the limit curve </> 2 <g> (j). But this is bent over in higher order in e 
as can be seen there since asymptotically the limit curve is reached. Generally one has to be 
careful with expansions like (jl,l|), (1.3) where higher orders are larger than previous orders for 



s » e/N. 

The O(N) and 0(<i)-scalar high-gradient operators leading t o (|1.3j ) have in fact already been 
investigated in an 1/iV-expansion by Vasil'ev and Stepanenko [[L2j ]. This remained unnoticed, 

3 As non-perturbative structures are always of interest in nontrivial field theories, our 1 /iV-results lead to the 
interesting conjecture that this limit structure holds to all orders in the iV-vector model. 
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probably due to the fact that powers of the auxiliary field a from (2.9) without any gradients 
turn out to interpolate between unstable high-gradient operators in 2 + e expansions and the 
unproblematic (c/> 2 ) s -operators in 4 — e expansions. This allowed the conclusion drawn from 
Eq. ( |4.6| ) that a stability problem for the nontrivial fixed point exists in first order in 1/N for 
s ~S> N below three dimensions. 

Of course one can argue that higher orders in the 1/iV-expansion will bend this behaviour over 
again, at least for some region of the <i-iV-plane. It is possible that the stability problem is only 
an artefact of the first order approximations and a proper treatment requires the investigation 
of the full perturbation series. This was just the situation for the limit point structure in the 
2 + e expansion discussed above. 



Considering the fact that other criticism of 2 + e or 1/iV-expansions |21, 22] also remained 
unresolved and the importance of 2 + e expansions in general, the present situation must be 
considered unsatisfactory. One possible direction for future work could be along the lines of 



Dasgupta and Lau [23, |24[. Their conclusion that topological excitations are relevant in the low 
temperature phase of the Heisenberg model could be linked to the high-gradient problem Q. 

Finally we pointed out in Sect. 4.2 that powers of the auxiliary field a lead to negative contact 
exponents ([4.1 1|) for 2 < d < 3. Negative contact exponents are very unusual in field theories 



and could be linked to multifractal scaling phenomena as argued in Ref. [19|. However, this 
interesting possibility relies on a stable nontrivial fixed point and hence a satisfactory solution 
of the above stability problem. 

This work was supported by Grant 95-01-00569a of the Russian Fond for Fundamental Research 
and by INTAS Grant 93-2492-ext and is carried out under the research program of the Inter- 
national Center for Fundamental Physics in Moscow. S.K.K. acknowledges hospitality of the 
Theory Group at St. Petersburg State University where part of this work was performed. 



Note added: 



After completion of this work we were made aware of Ref. p5[| where the role of topological 
defects in the d = N = 3 Heisenberg model has been re-examined using Monto Carlo simulations 
yielding different conclusions than Refs. [23, 24]. For recent work on topological defects in two 
dimensional systems see also Ref. [26|. 



Appendix A 



Here we show that under transformation ( |2.6D Eq. ( 2.18|) transforms into Eq. ( 2.21 ). It is 
obvious that it is sufficient to do this for the two particle case. For simplicity we consider the 
operator h\ 2 only 



1 pi— a 

da / d/3(l - a- f3) A ~ 2 ip(/3zi + (1 - a)z 2 , az 2 + (1 - &)z\) . 
o Jo 



(A.l) 



Then from the definition of the ^-function one has 



i>h(z 1 ,z 2 )= dt 1 dt 2 t 1 t 2 tl>h(ziti, z 2 t 2 )e 
Jo 



-ti-t 2 



(A.2) 
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After a change of variables (ti + *2 = T, ti = Ts) we get 

roc p\ pi rl—cz 

Mzi,z 2 ) = / dTT^e-T / d 7 / da / d/?( 7 (l - 7 )) A_1 

JO JO JO JO 

x (1 - a - (3) A - 2 iP(T((3 1Zl + (l - a )(l - 7)^2), T(a(l - 7)^2 + (1 - 0)7*i))- (A.3) 
We perform the following change of variables in the integral: 

7 = us - (1 - 7 )(1 - a) = (1 - a(l - 7) = t(l - v)\ (1 - /?h = (1 - *)(1 - f) (A.4) 

The last equation obviously follows from the three other equations. Then it is easy to derive the 
following relations: 

< 1/ < 1; 0<s + t<l (A. 5) 

and 

7(1 — rfdadfid'y = v(\ — v)dvdsdt , (A. 6) 

(1 _ a _ /?)A-2 (7(1 _ 7)) A-2 = (1 _ s _ t )A-2 (l/(1 _ i/)) A-2_ (AJ) 
Combining everything one immediately finds 

^h{zi,z 2 )= [ ds f dt(l-s-t) A - 2 4>(sz 1 + (l-s)z 2 ,tz 2 + (l-t)z 1 ) = h 12 4>(z 1 ,z 2 ). (A.8) 

JO JO 
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Figure captions 



Fig. 1. Cancellation of the self-energy diagram contribution of the cr-field. 

Fig. 2. Only diagram contributing in first order in 1/N for operator insertions from C sym . 

Fig. 3. The eigenvalues A obtained by diagonalizing ( |2.21| ) in the subspace of conformal operators 
with n = 3 fields plotted as a function of d = 2fi. The resp. critical exponents and 



anomalous dimensions can be obtained via Eq. (2.20). The number of gradients I is only 
given for the curves approaching the dashed limit curve 4> 2 (g> <j), the unlabelled curves 
approach other limit curves in agreement with our theorem. Notice that the normalization 
factor 2/r(d/2 — 2) has been chosen such that the values for d = 2 agree with the eigenvalues 
A2+ e in our definition ([2.281) , same for d = 4 and A4_ e using ( 2.29| ), 



Fig. 4. The function K(d) from fl4.6| ) plotted as a function of d. Observe the sign change in three 
dimensions separating stable from seemingly unstable behaviour. 
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